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Apart from omnidirectional, a solid elastic sphere is a multimode and multifrequency device 
for the detection of Gravitational Waves (GW). Motion sensing in a spherical GW detector thus 
naturally requires a multiple set of transducers attached to its surface at suitable locations. If these 
transducers are of the resonant type then their motion couples to that of the sphere, and the joint 
dynamics of the system has to be properly understood before reliable conclusions can be drawn 
from its readout. In this paper we address the problem of the coupled motion of a solid elastic 
sphere and a set of resonators attached to its surface in full theoretical rigor. A remarkably elegant 
and powerful scheme is seen to emerge from the general equations which shows with unprecedented 
precision how coupling takes place as a series function of the small "coupling constant" rj, the ratio 
of the resonators' average mass to the sphere's mass. We reassess in the new light the response of the 
^ , highly symmetric truncated icosahedron layout (the TJGA), and also present a new proposal (the 

_ O ■ pentagonal hexacontahedron, or PHC), which has less symmetry but requires only 5 rather than 6 

' transducers. We finally address the question of how the system characteristics are affected by slight 

departures from perfect spherical symmetry and identity of resonators, and find it to be quite robust 
against such small failures. In particular, we recover with fully satisfactory accuracy the reported 
experimental frequencies of the reduced scale prototype detector constructed and tested at LSU. 
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O ' I- INTRODUCTION 

■ 

I The idea of using a solid elastic sphere as a gravitational wave (GW) antenna is almost as old as that of using 
O . cylindrical bars: as far back as 1971 Forward published a paper in which he assessed some of the potentialities offered 
'' by a spherical solid for that purpose. It was however Weber's ongoing philosophy and practice of using bars which 
$_H . eventually prevailed and developed up to the present date, with the highly sophisticated and sensitive ultracryogenic 
systems currently in operation — see Q and for rather detailed reviews and bibliography. With few exceptions 
J> [Q;^)? spherical detectors fell into oblivion for years, but interest in them strongly re-emerged in the early 1990's, and 
an important number of research articles have been published since which address a wide variety of problems in GW 
spherical detector science. At the same time, international collaboration has intensified, and prospects for the actual 
I construction of large spherical GW observatories (in the range of ~100 tons) are being currently considered in several 
" " '" countries ^ , even in a variant hollow shape . 

A spherical antenna is obviously omnidirectional but, most important, it is also a natural multimode device, i.e., 
when suitably monitored, it can generate information on the GW amplitudes and incidence direction , a capability 
which possesses no other individual GW detector, whether resonant or interferometric. Furthermore, a spherical 
antenna could also reveal the eventual existence of monopole gravitational radiation, or set thresholds on it The 
theoretical explanation of these facts is to be found in the unique matching between the GW amplitude structure 
and that of the sphere oscillation eigenmodes: a general metric GW generates a tidal field of forces in an elastic 
body which is given in terms of the "electric" components Roioj (t) of the Riemann tensor at its centre of mass by the 
following formula |lC| ]: 

fGw(x,i) = Yl f^""Hx)g('™Ht) (1.1) 

1^0 and 2 
m — — 

where f ('""^ (x) are pure "tidal form factors", while g^^"^\t) are suitable linear combinations of the Riemann tensor 
components i?oiOj (t) which carry all the dynamical information on the G Ws monopole (I = 0) and quadrupole (/ = 2) 
amplitudes. 

On the other hand, a free elastic sphere has two families of oscillation eigenmodes, so called toroidal and spheroidal 
modes, and modes within either family group into ascending series of ?-pole harmonics, each of whose frequencies is 
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(2Z+l)-fold degenerate — see [0 for full details. It so happens that only monopole and/or quadrupole spheroidal 
modes can possibly be excited by an incoming met ric GW and their GW driven amplitudes are directly propor- 
tional to the wave amplitudes (jr^'™) (t) of equation (1.1). It is this very fact which makes of the spherical detector such 
a natural one for GW observations |1C| ]. In addition, a spherical antenna has a significantly higher absorption cross 
section than a cylinder of like fundamental frequency, and also presents good sensitivity at the second quadrupole 
harmonic 

But GW excitations are extremely weak | |T^ , and so a suitable readout system must be added to the sphere in 
order to monitor its motions and quantitatively assess their magnitude and physical significance. In cylindrical bars, 
current state of the art technology is based upon resonant transducers |p^-p^ . A resonant transducer consists in a 
small (compared to the cylinder) mechanical device possessing a resonance frequency accurately tuned to that of the 
cylinder. This frequency matching causes back-and-forth resonant energy transfer between the two bodies (bar and 
resonator), which results in turn in amplified oscillations of the smaller resonator. The energy amplification factor is 
-^resonator /-^bar, hence the amplitude amplification is the square root of this mass ratio. Pre-electronics mechanical 
amplification is highly desirable, given the exceedingly small magnitude of any expected GW signals arriving in the 
observatory. 

The philosophy of using resonators for motion sensing is directly transplantable to a spherical detector, but a multiple 
set rather than a single resonator is required if its potential capabilities as a multimode system are to be exploited 
to satisfaction. The practical feasibility of a multiple transducer readout system has been recently demonstrated 
experimentally with encouraging success by S. Merkowitz and W. Johnson at LSU, where a 740 kg prototype, milled 
in the shape of a truncated icosahedron, and endowed with a highly symmetric set of 6 resonators, was put to test 
[^-^. Their authors call this system TIG A, an acronym for Truncated Jcosahedron Gravitational Antenna. 

One of the pillars of the success of this prototype experiment has been M&J's ability to give an adequate theoretical 
interpretation of their experimental results. M&J's model of the coupled dynamics of sphere and resonators is based 
upon the hypothesis that only quadrupole excitations of the sphere's modes need to be considered in a GW detector, 
even if resonators are attached to the sphere's surface. This extrapolation of a result which, as we have just mentioned, 
does hold exactly for a free sphere, produces remarkably accurate predictions of the coupled system behaviour, too, 
even though interactions of the resonator set with the other, non-quadrupolc sphere's modes are neglected from the 
beginning. This is a strong indication that such neglected effects are second order for the accuracy of the experimental 
data. 

But what does "second order" precisely mean? The answer to such question requires the construction of a more 
elaborate model, which should be suitable to address in a systematic way any dynamical effects, and to quantitatively 
assess their real importance. The interest of a more sophisticated analysis is to understand and make clear the nature 
of a given approximation scheme, as well as to enable further refinement of it, if eventually required; its practical 
relevance is related to the reasonable expectation that future real spherical GW detectors will make use of extremely 
precise measurement techniques, likely to be rather demanding as regards accurate theoretical modeling of the system. 

The purpose and motivation of this paper is to present and develop such a more refined mathematical model. We 
address the the problem of the joint dynamics of a spherical elastic solid endowed with a set of radial resonators, 
with as few as possible unwarranted hypotheses, and with the objective to determine the system response to any 
interesting signals, whether GWs or calibration inputs, with unlimited mathematical precision. As we shall see, the 
model confirms and generalises that of Merkowitz and Johnson [ p^j22| in the sense of making precise its actual range 
of applicability. More specifically, we shall see that the solution to our general equations of motion can be written 
as a perturbative series expansion in ascending powers of the small coupling constant 77^/^ (j? = -MrosonatorZ-Msphcrc), 
whose lowest order terms exactly correspond to Merkowitz and Johnson's model. This is a key result, showing that 
the above alluded "second order" effects are precisely order effects, or that M&J's model is accurate up to relative 
errors of this order. 

Beyond this, though, a remarkably elegant, simple, and powerful algebraic scheme will be seen to emerge from the 
theory, which neatly displays the basic structure of the system for completely general resonator distributions over the 
sphere's surface. Based on the resulting equations we already advanced in references ||2^,Q a genuine transducer 
layout with 5 rather than 6 resonators, which we propose to call PHC (for pentagonal hexacontahedron, the shape 
of the underlying polyhedron), and which we shall also consider here in parallel with the more symmetric TIGA of 
Merkowitz and Johnson. 

The paper will be structured as follows. Section 2 is devoted to present the main hypotheses of the model and the 
general equations, whose reduction to a tractable set is described in detail in section 3. Sections 4 and 5 contain a 
number of general, signal-independent results which apply to an idealised system of perfect sphere and resonators; 
they constitute the foundations for the analysis of more realistic instances, where one (or more) of the idealised 
hypotheses partly fails. The system response to an incoming GW is then addressed in section 6, where a complete 
study of quadrupole and monopole radiation sensing is presented for arbitrary transducer layouts. Still within the 
section, we apply our general results to the important examples of TIGA and PHC. Section 7 addresses the problem 
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of the system response to a simple calibration signal, an impulsive hammer stroke, and in section 8 we assess the 
consequences of system deffects relative to the idealised perfection assumed that far; more specifically, we consider 
failures in spherical symmetry and identity of resonator properties, and tolerances in resonator location. In this section 
we also put to test our model's predictions by confronting them with the reported experimental data obtained in the 
TIG A prototype experiment [pl|, and sec that agreement between both (theory and experiment) is fully satisfactory 
in the given theoretical and experimental conditions. The paper closes in section 9 with a summary of conclusions. 

II. GENERAL EQUATIONS 

With minor improvements, we shall use the notation of references pO{ and [p3[|, some of which is now briefly 
recalled. We consider a solid sphere of mass M, radius R, (uniform) density g, and elastic Lame coefficients A 
and /i, endowed with a set of J resonators of masses Ma and resonance frequencies Qa (a = 1,. . . , J), respectively. We 
shall model the latter as point masses attached to one end of a linear spring, whose other end is rigidly linked to the 
sphere at locations Xa — see Figure |l|. The system degrees of freedom are given by the field of elastic displacements 
u(x, t) of the sphere plus the discrete set of resonator spring deformations Za(t); equations of motion need to be 
written down for them, of course, and this is our next concern in this section. 



a 



M 



q(t) 



Z(t) 



FIG. 1. Schematic diagramme of the coupling model between a solid sphere and a resonator. The notation is that in the text, 
but subindices have been dropped for clarity. The dashed-dotted arc line on the left indicates the position of the undeformed 
sphere's surface, and the solid arc its actual position. 

We shall assume that the resonators only move radially, and also that Classical Elasticity theory is sufficiently 
accurate for our purposesFI. In these circumstances we have [E3| 



= mV'u + (A + Ai) V(V-u) + f (x, t) (2.1a) 

Za{t)=-nl[Za{t)-Ua{t)]+^a{t) , a=l,...,J (2.1b) 

where ria = ^a/R is the outward pointing normal at the the a-th resonator's attachement point, and 

Uait) = na-u{Ka,t) , a=l,...,J (2.2) 

is the radial deform ation of the sphere's surface at x^. A dot (') is an abbreviation for time derivative. The term 
in square brackets in ( 2.1b ) is thus the spring deformation — q{t) in Figure |^. 



^ We clearly do not expect relativistic motions in extremely small displacements at typical frequencies in the range of 1 kHz. 



3 



f (x, t) in the rhs of (2.1a) contains the density of all non-internal forces acting on the sphere, which we expediently 
split into a component due the resonators' back action and an external action proper, which can be a GW signal, a 
calibration signal, etc. Thus 

f (x, t) — frcsonators ^) ~t~ fcxtcrnal(Xi (^■^) 



Finally, £,ait) in the rhs of ( |2.1b| ) is the force per unit mass (acceleration) acting on the a-th resonator due to 
external agents. 

Since we are making the hypothesis that the resonators are point masses the following holds: 

.7 

^resonators (^; ^) -Xa) [Za{t) - Ua{t)] Ua (2.4) 

a=l 

where J^^-* is the three dimensional Dirac density function. 

The external forces we shall be considering in this paper will be gravitational wave signals (of course!) and a simple 



calibration signal, a perpendicular hammer stroke. GW driving terms, we recall from (1.1), can be written 



fGw(x,i) = f("°)(x)5(°"Ht) + J2 f^""nx)5^^'")(t) (2.5) 

m=-2 

for a general metric wave — see ||l^ for explicit formulas and technical details. While the spatial coefficients f(''"^(x) 
are pure form factors associated to the tidal character of a GW excitation, it is the time dependent factors g^''™\t) 
which carry the specific information on the incoming GW. The purpose of a GW detector is to determine the latter 
coefficients on the basis of suitable measurements. 

If a GW sweeps the observatory then the resonators themselves will also be affected, of course. They will be driven, 
relative to the sphere's centre, by a tidal acceleration which, since they only move radially, is given by 

^GW(^) = c2 Ro^OJ{t) Xa,^na,J , = 1, . . . , J (2.6) 

where i?oiOj (0 are the "electric" components of the GW Riemann tensor at the centre of the sphere, and R is the 
sphere's radius so that, clearly, x^ — Rua- It may now be recalled from reference (lOj that 

c^i?o,o,(t)= E (2.7) 

1=0 and 2 
m— — 

where Efj^^ is a set of 6 (constant) symmetric matrices which verify^ 

Efj^Krnj^Yi„,{n) , 1 = 0,2, m^-l,...,l (2.8) 

Finally thus: 

^^'^{t)^R E yim{na)9^'"'\t) , a=l,...,J (2.9) 

;=0 and 2 
m— — 

We shall also be eventually considering in this paper the response of the system to a particular calibration signal, 
consisting in a hammer stroke with intensity fp, delivered perpendicularly to the sphere's surface at point Xq: 



^stroke 



(x,t)=fo(5(3)(x-xo)<5(i) (2.10) 



which we have modeled as an impulsive force in both space and time variables. Unlike GW tides, a hammer stroke 
will be applied on the sphere's surface, so it may have no direct effect on the resonators. In other words, 

e"''W=0, a=l,...,J (2.11) 



Y;m(n) are spherical harmonics 
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Our fundamental equations thus read 



92 u 



J 

E 

b=l 



Mb^l 5(3) (X - Xb) [Zb{t) - Ub{t)] rib + foxtornal(x, t) 



Zait) = [Za{t) - Ua{t)] + Ut) , a=l,...,J 



(2.12a) 



(2.12b) 



where fcxtornai(x, t) wiU be given by either ( ^.Sj ) or ( 2.10| ), as the case may be. Likewise, <^a(i) wiU be given by ( |2.9| ) 
or ( ^.11 ) , respectively. The remainder of this paper will be concerned with finding solutions to the system of coupled 
differential equations ( 2.12 ), and with their meaning and consequences. 



III. GREEN FUNCTION FORMALISM 



An elegant and powerful method to solve equations (2.12) is the Gree n fu nction formalism. It so happens that, 



in all instances of our concern here, the force density f(x, t) of equation (2.3) is of the separable type, i.e., it can be 
written as a sum of pro duc ts of a function of the space variables x times a function of the time variable t. Direct 
inspection of equations (2.4)-(2.10) readily shows that this is always the case. We thus have, generically. 



f(x,t)=^f(")(x)gW(t) 



(3.1) 



where a is a suitabl e labe l. We recall from reference ||T^ that, in such circumstances, a formal solution can be written 
down for equation ( 2.1a ) in terms of a Green function integral, whereby the following orthogonal series expansion 
obtains: 



u(x,t) =5]^ c.^V^"^u^(x)g(,")(t) 



a N 



where 



f(a) 

Jn 



/ u^(x).f(")(x)d3 



g^^\t)^ / g(")(t') smu;^{t^t')dt' 



(3.2) 



(3.3a) 



(3.3b) 



Here, lun and ujv(x) are the eigenfrequencies and associated normalised wavefunctions of the free sphere — see 
again ||l^ for a comprehensive characterisation. Also, N is an abbreviation for a multiple index {nlm}. We quote the 
result of a few explicit calculations which will be useful later on: 



/\cij q2 * ^ 1 7 

resonators, Af — ^ 'a "^a'^Jvl^aj , O— i,...,J 

fGW,N = ^11' <^™™' ' ^ = {nlm} , /' = 0, 2 , 

/stroke, W = ^ fo ' (xq ) 



Aa) 



A' ...A' 



"resonators , 



where the coefficients a„; in ( 3.4b ) are overlapping integrals of f'^'™)(x) across the volume of the sphere 

r{t)= f [Za{t')~Uait)]smLUNit-t')dt' , a=l,...,J 

9^w\it)^ I g^'^^^t') sinujN {t-t')dt' 
Jo 

ffstroke,Ar(i) = SlU UJjyt 



(3.4a) 

(3.4b) 
(3.4c) 
, and 

(3.5a) 

(3.5b) 
(3.5c) 
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These expressions can now be substituted into equation (^^) to obtain the following equivalent set of equations of 
motion: 



U(x,t) = "^N \ Y^^l [nb-U^(xb)] 5rcLator.,A'(i) + /ixtornal.JV slxternal.AT (*) 



N kfc=l 
12 



M 



Za{t) = ~K [Za{t) - Ua{t)] + , = 1, . . . , J 

We now specify x = Xa in ( |3.6a ) and multiply both sides of the equation by n^, thus finding 

Ua{t) = ^r^-^^HO + E / ^-"(^ - ^'^b{t') - Ub{t')] dt' 

Za{t) = -nl [Za(t)-Ua{t)]+^a{t) , a ^ I, . . . , J 

where uexternaij-^^ = n„ • u'='''°™^i (x^ , i), and 



(3.6a) 
(3.6b) 

(3.7a) 
(3.7b) 

(3.8) 



a N 



is the sphere's response to an external force in the absence of resonators . The form of (3.8) is given in reference 
pO| for a generic metric GW signal and for a hammer stroke. The kernel matrix Kabit) in ( 3.7b| ) is the following 
weighted sum of diadic products of wavefunctions: 



E 

N 



[n6-u^(x,,)] [nQ-UAr(xa)] sinwjvi 



Finally, we have defined the mass ratios of the resonators to the entire sphere 



lb 



Mb 
M 



.J 



(3.9) 



(3.10) 



which will be s mall parameters in a real device. 

Equations (3.7) are a set of integro-differential equations for the radial deformations of the sphere, Ua{t), and those 
of the resonators, Za{t). If we solve them t hen w e obtain at on ce the complete solution to our general problem by 
direct substitution of these quantities into (3.5a), then in ( |3.6a| ). But before going into the technical details of the 
solution pro cess l et us briefly pause for a qualitative inspection. 

Equation ( 3.7a ) shows that the sphere's deformations Ua{t) are made up of two contributions: one due to the action 
of external agents (GWs or other), contained in w|j'^*°''"'''(i), and another one due to coupling to the resonators. The 
latter is commanded by the small parameters rjb, and correlates to all of the sphere's spheroidal cigcnmodcs through 
the kernel matrix Kab{t). This has consequences for GW detectors, for even though GWs only couple to quadrupole 
and monopole]^ spheroidal modes of the free sphere [^,^ 11 1, attachement of resonators causes, as we see, a well 
defined amount of energy to be transferred from these into other modes of the antenna, and conversely, these modes 
back-act on the former. As we shall shortly prove, such effects can be minimised by suitable tuning of the resonators' 
frequencies, but outright neglection of them results in inaccurate conclusions about the system dynamics. 



A. Laplace transform domain equations 



We now take up the problem of solving equations (p.7|). Equation ( 3.7a ) is an integral equation belonging in the 
general category of Volterra equations |29[ | , but a series solution to it in ascending powers of rjb by iterative substitution 
of Ub{t) into the kernel integ ral is not viable due to the dynamical contribution of Zb{t), which is in turn governed 
by the differential eq uati on ( 3.7b ). A better suited method to address this integro-differential system is to Laplace- 
transform equations (3.7). We denote the Laplace transform of a generic function of time f{t) with a caret (") on its 
symbol, e.g.. 



Monopole modes only exist in scalar-tensor theories of gravity, such as e.g. 
belong in this category. 



Brans-Dicke 



General Relativity does not 
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Jo 

and make the assumption that the system is at rest before an instant of time, t = Q, say, or 

u(x,0)-u(x,0) = z,(0) = i,(0) = 



Equations (3.7) are then recast in the equivalent form 



(3.11) 



(3.12) 



na{s) = ^r^^'^s) - ^ ry, kab{s) [h{s) - Ub{s) 



6=1 



Za{s) = [Zais) - Ma(s)] + ^a(s) , fl = 1, . . . , J 



(3.13a) 



(3.13b) 



for which use has been made of the convolution theorem for Laplace transformsQ. A further simplification is 
accomplished if we consider that we shall in practice be only concerned with the measurable quantities 



qa{t) = Za{t) ~ Ua{t) , a=l,...,J 



(3.14) 



representing the resonators' actual elastic deformations — cf. Figure It is readily seen that these verify the 
following: 



E 

b=l 



Sab + Vb ■ 



Kabis) 



qb{s) 



"■external 



where 



and 



--, external 



a L 



E 

N 



/ixternal, AT ["a ' (Xa ) 



(^) 



Us) 



a = 1, . . . , J 



5external(*) 



Kabis) 



N 



[nb-U^(xb)] [na-UAr(Xa) 



(3.15) 



(3.16) 



(3.17) 



which ensu e dire ctly from (3.8), (3.9) and the definition ( [3.11 ). 

Equations (3.15) constitute a significant simplification of the original problem, as they are a set of just J algebraic 
rather than integral or differential equations. We must solve them for the unkno wns Qajs), then perform inverse 
Laplace transforms to revert to qa{t)- But we note that Equations (3.15) and ( 3.16| ) indicate that the Laplace 
transform functions we shall be concerned with are of the rational class, i.e., they are quotients of polynomials in s. 
This greatly facilitates the latter step (Laplace transform inversion) , as it can in this case be obtained by the calculus 
of residues . 

Determination of the poles of (jai-s) is thus required in th e firs t place. Clearly, poles correspond to those values 
of s for which the matrix in square brackets in the Ihs of ( |3.15| ) is sm^it/ar^ or, equivalently, to the zeroes of its 
determinant: 



A(s) = det 



6ab + Vb ■ 



Kabis) 







Poles 



(3.18) 



As is well known, the imaginary parts of the poles are the system characteristic frequencies, or resonances 
residues at such poles determine the specific weight of the respectively associated modes in the system response to a 
given external agent. 

The procedure's general guidelines to solve the problem are thus clear-cut. The details of its actual implementation 
are however not obvious, so we now come to them. We begin with the simplest case, i.e., a perfect sphere with 
perfectly tuned identical resonators. 



* This theorem states, it is recalled, that the Laplace transform of the convolution product of two functions is the arithmetic 
product of their respective Laplace transforms. 
^ I.e., non-invertible. 
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IV. PERFECT SPHERE, IDENTICAL RESONATORS, IDEAL TUNING 



We first consider this highly ideahsed situation as a suitable starting point to address more realistic practical 
instances, which we naturally do not expect to depart excessively from that ideality, anyway. Moreover this simpler 
analysis will be very interesting conceptually, as all the major (coarse grain) characteristics of the system will emerge 
out of it. 

In a perfect sphere the frequencies of the eigenmodes can be classified as Z-pole series of ascending harmonics, each 
frequency in a given series being (2Z+l)-fold degenerate. GWs exclusively couple to spheroidal eigenmodes and 
our resonators are assumed to be sensitive to radial deformations of the sphere, so only spheroidal frequencies will 
concern us here. Following the notation of reference we denote by uini the n-th harmonic of the Z-th /-pole series, 
whose associated 2l-\-l eigenf unctions have the following radial projections: 



n-u„i„(x) ^ Ani{r)Yijn{e,ip) 



(4.1) 



where Yim{0,^) are spherical harmonics and Anijr) are given in |T^. Degeneracy of the eigenfrequencies ujni 
enables direct summation over the degeneracy index m in (3.17), so that 



^ \A,AR)\' - 



- Piii^a-rib) = ^ 



(nl) 

,2 ^ab 
'■'nl 



(4.2) 



for a perfectly symmetric sphere. Here, use has been made of the summation formula (Al) — see Appendix below. 
Our next assumption is that all the resonators are identical, or 



Vi = 



= nj = n 



(4.3) 



The fundamental idea behind using resonators is to have them tuned to one of the frequencies of the sphere's 
spectrum, so we now make our third hypothesis: 



(4.4) 



In a GW detector it will only make sense to choose — or Iq = 2, of course, and have no refer to the first or perhaps 
second harmonic We keep the generic expression (4.4) for the time being to encompass all the possibilities within 
a single formalism, and to make room for calibration signals, too. 

Based on the above thre e hypotheses — i.e., perfect spherical symmetry plus equations ( [4.3| ) and (4.4) — , we can 
now rewrite equation (3.17): 



E 

b=l 



7ll 



(nl) 
ab 



'■.external 



f]2 



is) 



Us) 



Equation ( |3.18 ) for the system resonances can also be recast in a more convenient form: 



(f^-c^„„,J (4.5) 



A(s) = det 



5ab + ?7 



(s2 +f]2)2 



("o'o) 



E 



(nl) 



{s^ + n^){s^ + Loi,) 



= 



(4.6) 



To find an analytic expression for the roots of (|4.6D is an impossible task, so we resort to a perturbative expansion 
in the small "coupling constant" rj. As we shall now see, much light is shed onto the physics of the problem as we 
proceed with that expansion. 



V. FREQUENCY SPECTRUM 

Let us formally state the essential fact that the resonators are much less masssive than the sphere, which is the 
basis for all the considerations from now on in this paper: 

77 <1 (5.1) 



It is clear from the structure of equation (4.6) that vanishing of A(s) will require s to be such that the denominators 
in the fractions within square brackets there be proportional to r?; otherwise it will not be possible to have A(s) — 
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for arbitrarily small values of r], since Sab is of course a regular matrix. But this points to a sharp distinction between 
roots which are close to s = zLi^l = zLiujnoioj ^-nd roots which are close to s = zLiUni {nl ^ nolo). The former's degree of 
proximity to ±if2 is of order while the latter's proximity to iiujni is of order 77, instead. We thus correspondingly 
distinguish the following two categories of roots: 

sl = -n^ {l + Xiv'^^ + XiV+--) {^ = ^noio) (5.2a) 
4i = -^li (1 + t'^ V + bi"'^ v' + --) {nl + nolo) (5.2b) 

The coefficients xLj Xiv and 6^"''', can be calculated recursively, starting form the first. As we now show, 

these two categories of roots present qualitatively different characteristics. We consider them separately. 



A. Roots near Q 



Upon substitution of (5.2a) into (4.6) it is readily seen that xj. is a solution to the algebraic equation 



det 



'^ab - -2 Xab 



= 



(5.3) 



So xi are the eigenvalues of the matrix xi&°'''\ which happens to be non-negative definite — see Appendix. There 

2 

are of course J of these, some of which may be repeated {degenerate) , or null. This means that we get J pairs of roots 
of the type snin ( |5.2a ) lying on the imaginary axis of the complex s-plane symmetrically (to order 7/^/^) around ztifl 
— see Figure |2[ The system resonant frequencies can thus be represented by the pairs 



1 ± 



2/ + 1 

An 



\An„UR)\ CaV'^']+0{7l) , a=l,...,J 



(5.4) 



where the eigenvalues of the matrix Pig{na-vib), and 0{r]) stands for the higher order terms in (5.2a) 



>< s 



iO. 



nl 



X s, 



>< s, 



s-plane 



FIG. 2. Qualitative situation of the roots in the complex s-plane. Dots correspond to the free sphere's frequencies, while 
crosses correspond to coupled system frequencies. Roots near the tuning frequency fl come in pairs (we draw just one, 
corresponding to one single resonator, to avoid unnecessary complication at this stage), while roots near other (high er) fr ee 
sphere frequencies are downshifted rather than split, though the effect is altogether weaker in this case — see equation (5.2b). 



Equation ( |3.4[ ) very much reminds us of a familiar result in cylindrical bar theory fl^ : attachement of a resonator 
to the bar's end face causes the latter's characteristic resonance frequency to split up into a symmetric pair around 
the original value, the amount of relative shift being proportional to 77^/^. As we now see, in a spherical body the 
number of frequency pairs equals the number of resonators, it is of the same order of magnitude (77^/^), and is precisely 
controlled by the geometry dependent eigenvalues C.a- 
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It must however be cautioned that ( |5.4| ) is not an exact result, but the d omin ant approximation in an ascending 
series in powers of ry^^^. If better accu racy is needed then higher terms in ( 5.2a ) must be calculated. For example, 
the next order correction to Sq in ( ^.2a ) is 



_ 1 



A.aa 



1,. 



(5.5) 



corresponding to the a-th root of the previous order, equation ( ^.3[ ). In this formula we have Xaa \ which stands for 
the a-th diagonal element of the matrix xib''' in the (abstract) system of axes in which x^'^"'"'' is diagonal. 



The second order corrections ( p.5[ ) are seen to be equal for the two members of the frequency pairs (5.4), as 
they only depend on Xi ■ They thus result in rigid shifts for both of them. More substantial, the actual value of 

2 

these corrections involves the whole spectrum of spheroidal eigenfrequencies and wavefunctions. This should not be 
considered surprising on general groundsj^ but to find this result has required us to model the system dynamics 
in a more elaborate conceptual structure than has been considered so far in the literature pO| , ^ , ^ . The physical 
interpretation of ( |5.5| ) is that the attachement of resonators causes a certain level of intercommunication between 
different vibration eigenmodes of the sphere, even if the resonators are accurately tuned to one of the spectral 
frequencies. This "cross talk" between modes has c onse quences for non-tuned modes, as we shall see in the next 
subsection, as well as for tuned modes, as we see in (5.5). The practical relevance of these higher order corrections 
will naturally be dictated by the precision of measurments in a real GW antenna. As of this date, the experimental 
data available on such devices as described here is to our knowledge limited to the TIGA prototype data |^ , but the 
experimental conditions reported [2^ are not sufhciently fine tuned to guarantee that corrections of order rj to fit the 
data are really meaningful — see section 8 below. This is surely the reason why errors in the formulas given by other 
authors for system frequencies Bq] have passed unnoticed so far. 



B. Roots near cuni 7^ ^ 



We now come to the other roots, equation (5.2b). If the latter is substituted into (^■(\ ) then we readily see that 
6^"''' satisfies the algebraic equation 



dot 



^"i ,inl) 



hr Sab 



(nl) 
Xab 



= 



(5.6) 



which is again an eigenvalue equation, having J solutions. Since the eigenvalues of the matrix xih'^ either positive 

or null (see Appendix) it follows that 6^"'^ either has the same sign as (£7^ — cj^^) or is zero, respectively. Thus for 
example, if Q is chosen equal to the lowest frequency of the spectrum (first quadrupole harmonic, UJ12) then the 
system frequencies close to tOni will be slightly downshifted relative to their original value coni. This is schematically 
represented in Figure ^. 

We see again in equation ( ^ ) that the presence of resonators does indeed affect the whole spectrum of the free 
sphere, though with an intensity of coupling which differs between tuned and non-tuned modes by factors of order 
7^1/2 shall not go into any more depth in the analysis of these fine structure details in this paper. 



VI. SYSTEM RESPONSE TO A GRAVITATIONAL WAVE 

Our next concern is the actual system response when it is acted upon by an incoming GW, i.e., which are the 
amplitudes of the excited modes, whose frequencies we have just estimated, and how do they relate to the GW 
amplitudes g^^"^\t). To this end we must find the inverse Laplace transform of qa{s)'. 



qa{s) 



E 

6=1 



. \p ^^S^ (nl) 

ab+V 2^ (s2 + n^){s^ + LU„i)^ 



-,GW 



S2 -f 



(^) 



Us) 



,J (6.1) 



Second order perturbative corrections to eigenvalues often show this feature - 
Mechanics llsl. 



-recall e.g. energy level corrections in Quantum 
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as follows from (|4.5|), where ■u^'^*''™'*'(s) has been substituted by ?ip^(s), which is in its turn given by 



=0 and 2 \n=l / 



(6.2) 



according to ( |3.4b| ), ( |3.5b| ) and (|^. Clearly, from (|2j) we also have 



a = 1, 



.,J 



(6.3) 



Z=0 and 2 



In order to ease the notation we can combine (3.1) and (3.2) into a more compact form: 

Us)= J2 <^("")(s)5('™)(s) , a=l,...,J 



(6.4) 



I = and 2 
m — —I, ...A 



with, obviously. 



$i'™)(.)^-^ 
f)=i 



'^E 



(s2+r!2)(s2+C^„i) 



2 -^ah 



-4 + E 



flni ^ni (-R) 



We thus have, by the convolution theorem again. 



1 Jo 

Im 



a — 1, 



YlrniUb) (6.5) 



(6.6) 



where ^a'^\t) is the inverse Laplace transform of ( [6.51 ). Despite its complexity, <&a™'*(s) is a rational function 
of s, and therefore its inverse Laplace transform can be calculated, as already advanced in section 3, by the residue 
theorem through the formula pl[| 



$(^'™)(t) = 27ri ^ jresidues of (s) e"* at its poles in complex s-plane| 



(6.7) 



The poles in this formula happen to be exclusively those of equations (5.2), and there are no poles either at = —fl^ 
or at = — w^;, for divergences compensate each other in numerators and denominators in (6.5) at those points. Our 
next step is thus to calculate the corresponding residues. This is a relatively simple, but considerably laborious task, 
the details of which will be omitted here. We just quote the results and directly move on to their discussion. 



One can see that ^a"^\t) has the following general structure: 



$i""Ht) oc 77-1/2 J2 (sinc^e+t - sinc^e-t) Suo + 0(0) 

Cc#o 



(6.8) 



where ujc± are given in ( p.4\ j, and the sum above excludes terms associated to null eigenvalues ^c- The term 0(0) 
stands for higher order corrections, which in this case happen to be of order zero in 77. But let us look at the important 
qualitative consequences of (|6.8[). 



i) The resonators' motions occur with a mechanical amplification factor of 77^^/2 relative to the driving GW 
amplitudes (t), provided the resonator frequency Q is tuned to a monopole or quadrupole sphere frequency 
— see the factor Su„ in ( |6.8| ). If other frequencies are chosen for tuning then there is no coupling to GWs to this 
order, but weaker by factors of 77^/2^ ^t least. 

i) The spectral composition of these motions is dominated by the symmetric frequency pairs Wci into which the 
tuned frequency (-Un„ia = ^ splits up as a consequence of the resonators' presence. The maximum number of 
such frequency pairs is (2Zo+l), e7;e7i if the number of resonators, J, is larger than this figure. This is because 
(2Z0+I) is the maximum number of non-null eigenvalues — cf. Appendix. 
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iii) Higher order terms, sym bolis ed by 0(0) in (6.8), encompass couplings to all the non-tuned modes, whose 
frequencies are given by ( 5.2b| ), as well as corrections of the same order of magnitude in the tuned mode. The 



latter include in particular modes corresponding to a null eigenvalue ^c, which are guaranteed to be present 
whenever J > 2lo+l. As we shall shortly see, an example of this is provided by the TIGA layout [pT| . 

Higher order corrections are considerably more difficult to address than first order, and we shall not attempt a 
closer approach to them in this paper. So let us continue with the investigation of the system response to lowest 
order. Given the above discussion, it is expedient to recast (^^) in the form 

ga(s)=ry-i/2 ^ A('™)(s;17)g('™)(s)+O(0) , a=l,...,J (6.9) 

Lrn 

where D, will be tuned to either a monopole (fl — ujno) or a quadrupole (ri = a;„2) frequency of the free sphere's 
spheroidal spectrum. We consider the two possibilities successively. 



A. Monopole gravitational radiation sensing 



General Relativity, as is well known, forbids monopole GW radiation. More general metric theories, e.g. Brans- 
Dicke |2^, do however predict this kind of radiation. It appears that a spherical antenna is potentially sensitive to 
monopole waves, so it can serve the purpose of thresholding, or eventually detecting them. This clearly requires that 
the resonator set be tuned to a monopole harmonic of the sphere, i.e., 



= w„o , ('o = 0) 

where n tags the chosen harmonic — most likely the first (n = 1) in a thinkable device. 
Since Po{z) = 1 (for all z) the eigenvalues of Po(na-nb) are, clearly, 

cf^j, ci^ .■■=cj = o 

for any resonator distribution. The tuned mode frequency thus splits into a single strongly coupled pair: 



(6.10) 



(6.11) 



0(7?) 



^ — w„o 



(6.12) 



The A-matrix of equation (f^M is easily seen to be in this case 



Ai'™)(.;..„o) = (-l)^^^ ' 



whence the system response is 



-1/2 



i-iy 



1 r 



5(°°)(s) + O(0) 



(6.13) 



(6.14) 



regardless of resonator positions. The overlap coefficient a„o is calculated by means of formulas given in ||l0|] , and has 
dimensions of length. By way of example, aio/i?— 0. 214, a nd a2a/R— —0.038 for the first two harmonics. 

A few interesting facts are displayed by equation ( |6.14 ). First, as we have already stressed, it is seen that if the 
resonators are tuned to a monopole detector frequency then only monopole wave amplitudes couple strongly to the 
system, even if quadrupole radiation amplitudes are significantly high at the observation frequencies lu±. Also, the 
amplitudes cjais) are equal for all a, as corresponds to the spherical symmetry of monopole sphere's oscillations, and 
are proportional to J"-'^/^, a factor we should indeed expect as an indication that GW energy is evenly distributed 
amongst all the resonato rs. A single transducer suffices to experimentally determine the only monopole GW amplitude 
^(oo)|-g-j^ of course, but ( |6.14 ) provides the system response if more than one sensor is mounted on the antenna for 
whatever reasons. 
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B. Quadrupole gravitational radiation sensing 



We now consider the more interesting case of quadrupole motion sensing. We thus take 

n = w„2 , (^0 = 2) 



(6.15) 



where n labels the chosen harmonic — most likely the first (n= 1) or the second {n = 2) in a practical system. The 
evaluation of the A-matrix is now considerably more involved, yet a remarkably elegant form is found for it: 



6=1 



(c) (c)* 



Cc 



Y2minb)Sl2 (6.16) 



where vi'^^ is the c-th normalised eigenvector of P2(na-nb), associated to the non-null eigenvalue Cc- Let us stress 
once more that equation (6.16) explicitly shows that at most 5 pairs of modes, of frequencies uJc±, couple strongly 
to quadrupole GW amplitudes, no matter how many resonators in excess of 5 are mounted on the sphere. The tidal 
overlap coefhcients a2„ can also be calculated, and give for the first two harmonics [E3[ 



Ql2 

R 



= 0.328 



022 

R 



0.106 



(6.17) 



The system response is thus 



9q(s) = -q 



b=i yca^o 



(s'+u^i^y'^is' + ^n 



(c) (c)* 



X V y2™(nb)g(2'")(s)+O(0) 



E 

m=-2 



(6.18) 



Equation ( |6.18 ) is completely general, i.e., it is valid for any resonator configuration over the sphere's surface, and 
for any number of resonators. It describes precisely how all 5 GW amplitudes g^^'"-'(s) interact with all 5 strongly 
coupled system modes; like before, only quadrupole wave amplitudes are seen in the detector (to leading order) when 
^^ = '^n2, even if the inco ming wave carries significant monopole energy at the frequencies ujc±. 

The consequences of (3.18) are best seen in practical examples, so we now come to a more detailed consideration 
of two specific resonator distributions. 



C. The TIG A configuration 

A highly symmetric resonator layout has been proposed and experimentally studied by Merkowitz and Johnson at 
LSU ||l3-Gl|, which consists in a set of six transducers attached to the six non-parallel pentagonal faces of a truncated 
icosahedron, as shown schematically in Figure ^. 

It can be immediately verified that the five non-null eigenvalues of P2{na-nb) are all equal for this distribution, and 
there is a null sixth, too: 

^2 = • • • = C2' - ^ , = (TJGA) (6.19) 

This means that all five pairs of strongly coupled frequencies collapse into a single, five-fold degenerate pair: 

Lul = tol2 (1 ± \An2{R)\ V'^^^ + 0{v) , a = 1, . . . , 6 (6.20) 

Degeneracy in this layout is a consequence of its symmetry relative to the quadrupole structure of the GW driving 
force and the tuned detector modes. It is also not difficult to see that the TIGA is the minimal configuration with 
so much degeneracy, as there are no 5 resonator configurations with equivalent symmetry. There are however other 
non-minimal sets with the same degree of degeneracy — for example 10 resonators on the ten non-parallel faces of a 
regular icosahedron, etc., see e.g. [p3l for further analysis of this and other possibilities. 
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FIG. 3. Schematic view of the TIGA: a truncated icosahedron shape with black dots at the resonators' posotions. 



The normalised eigenvectors associated to the five non- vanishing eigenvalues (3.19), which are the only relevant 
ones now, are seen to be 



,2 , a= 1,, 



,6 



(TIGA) 



The system response is thus given by 



9q(s) = V 



- ^-1/2 



1 r 



3 2 



'] ^ r2™(n.)g(2'")(s) + O(0) , a = l 



(6.21) 



(6.22) 



m=-2 



This highly symmetric and remarkably simple formula was obtained for the first time by Merkowitz and Johnson 
1^,^. Its scope and range of validity, as well as its uniqueness, are now more firmly established in the light of the 
present, more elaborat e ana lysis. 

Based on equations ( 6.22| ) Merkowitz and Johnson define what they call mode channels: these are linear combina- 
tions of the resonators' readouts which directly yield the GW amplitudes (j^'^ "^Hs) at the frequencies uj±, and they 
are easily obtained thanks to the orthonormality property of the eigenvectors ( 6.2l| ). They are 



y 



(m)f 



(^)^^„(™)*^^(,) = ^-i/2„^^ 1 r(,2^^2) '^^s' + ^iy']gi^-^)^,) + o{0), m = -2,...,2 (6.23) 



There naturally are 5 rather than 6 mode channels — as there are 5 quadrupole GW amplitudes — , and it is easy 
to implement an algorithm to calculate them on line from raw detector data. So signal and direction deconvolution 
methods [sl JlOlJs^ can be directly applied to the mode channel set very adva ntage ously. This should be considered one 
more attractive property of the TIGA layout since, as we see in equation (S.18), arbitrary resonator configurations 
do not generally permit the construction of such mode channels so efficiently. 



D. The PHC configuration 



As we have just seen the TI layout is highly symmetric, and is the minimal set with maximum degeneracy. To 
accomplish this, however, 6 rather than 5 resonators are required on the sphere's surface. Since there are just 5 
quadrupole GW amplitudes one may wonder whether there are alternative layouts with only 5 resonators. Equation 
( |6.18 ) is completely general, so it can be searched for an answer to this question. In reference we made a specific 
proposal, which we now describe in more detail. 

In pursuing a search for 5 resonator sets we found that distributions having a sphere diameter as an axes of 
pentagonal symmetri^ exhibit a rather appealing structure. More specifically, let the resonators be located at the 
spherical positions 



By this we mean resonators are placed along a parallel of the sphere every 72°. 
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27r 

9a ^ a (alia), (^a = (a"l)^, a = 1, . 

5 



The eigenvalues and eigenvectors of P2(na-nfc) arc then 



Co^^l (3 cos^a-l)^ 



Ci = C-i = ¥ sin"' a cofp- a 



CI = e 



sm a 



= y^C™'>"2rn(na) , m = -2, . . . , 2 , a = l,...,5 
so the A-matrix is also considerably simple in structure in this case: 

Ai'™)(s;c^„2) = 
where we have used the obvious notation 



^1 /--I 1 

J. 0'n2 <5„i 2 



(' 2 I 2 



'-(s^+c^i 



>2m(na)<5z2 (PHC) 



(6.24) 

(6.25a) 
(6.25b) 

(6.26) 



■^"2 I 1± |A.2(i?)| Cm?? 



1/2 



0(7?) 



m = -2,...,2 



(6.27) 



As we see from these formulas, the five expected pairs of frequencies actually reduce to three, so pentagonal 
distributions keep a certain degree of degeneracy, too. The most important distinguishing characteristic of the general 
pentagonal layout is best displayed by the explicit system response: 



qa{s) = -V 



^ _„-l/2 ^ /4ir 



1 

1 

1 
2^ 



(6.28) 



This equation indicates that different wave amplitudes selectively couple to different detector frequencies. This 
should be considered a very remarkable fact, for it thence follows that simple inspection of the system readout 
specirMTTT^ immediately reveals whether a given wave amplitude g^™'{s) is present in the incoming signal or not. 

Pentagonal configurations also admit m ode ch annels, which are easily constructed from ( 6.28| ) thanks to the or- 
thonormality property of the eigenvectors ( |6.25b| ): 



^,W(.) ^ -i"^*9<^(-) = a„2 ^ (.s^ + u:^)-' - (.^ + g('"ys) + 0(0) 



(6.29) 



These are almost identical to the TIGA mode channels (6.23), the only difference being that each mode channel 
comes now at a single specific frequency pair uj,n± . It thus appears that a pentagonal transducer configuration enables 
signal observations over a somewhat richer frequency band than does the TIGA. 

Based on these facts one may next ask which is a suita ble tr ansducer distribution with an axis of pentagonal 
symmetry. In Figure ji] we give a plot of the eigenvalues ( 6.25a ) as a function of a, the angular distance of the 
resonator set from the symmetry axis. Several criteria may be adopted to select a specific choice in view of this 
graph. An interesting one was proposed by us in reference |23|| with the following argument. If for ease of mounting, 
stability, etc., it is desirable to have the detector milled into a close-to-spherical polyhedric shape^ then polyhedra 
with axes of pentagonal symmetry must be searched. The number of quasi regular convex polyhedra is of course finite 
— there actually are only 18 of them |3^,^^ — , and we found a particularly appealing one in the so called pentagonal 
hexacontahedron (PHC), which we see in Figure |5[ left. This is a 60 face polyhedron, whose faces are the identical 
irregular pentagons on the right of the Figure. The PHC admits an inscribed sphere which is tangent to each face 



In a noiseless system, of course 
^ This is the philosophy suggested and experimentally implemented by Merkowitz and Johnson at LSU. 
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FIG. 4. The three distinct eigenvalues (^rn ("i = 0,1,2) as funct ions o f the distance of the resonator parallel's co-latitude a 
relative to the axis of symmetry of the distribution, cf. equation (6.25a). 



at the central point marked in the Figure. It is, clearly, to this point that a resonator should be attached so as to 
simulate an as perfect as possible spherical distribution. 

The PHC is considerably spherical: the ratio of its volume to that of the inscribed sphere is 1.057, which quite 
favourably compares to the value of 1.153 for the ratio of the circumscribed sphere to the TI volume. If we now 
request that the frequency pairs be as evenly spaced as possible, compatible with the PHC face orientations, 
then we must choose a = 67.617°, whence 



ujo± = ijjyi ( 1 ± 0.5756 r\ 



1/2 



2 (l± 0.8787 ?7 



1/2 



^2± = ^\ 



2 (l ± 1.0668 (6.30) 



for instance for 17 = ^12, the first quadrupole harmonic. In Figure ^ we display this frequency spectrum together 
with the multiply degenerate TIG A for comparison. The frequency span of both distributions is naturally comparable, 
yet the PHC is slightly broader. 

The criterion leading to the PHC proposal is of course not unique, and alternatives can be considered. For example, 
if the 5 faces of a regular icosahedron are selected for sensor mounting (a = 63.45°) then a four-fold degenerate pair 
plus a single non-degenerate pair is obtained; if the resonator parallel is 50° or 22.6° away from the "north pole" then 
the three frequencies wo+, 1^1+, and 0^2+ are equally spaced; etc. The number of choices is virtually infinite if the 
sphere is not milled into a polyhedric shape. 

Let us finally recall that the complete PHC proposal was made with the idea of building an as complete as 
possible spherical GW antenna, which amounts to making it sensitive at the first two quadrupole frequencies and 
at the first monopole one. This would take advantage of the good sphere cross section at the second quadrupole 
harmonic p^ , and would enable measuring (or thresholding) eventual monopole GW radiation. Now, the system 

pattern matrix Aa™-*(s;f2) has identical structure for all the harmonics of a given / series — see ( 6.13[ ) and (3.16) — , 
and so too identical criteria for resonator layout design apply to either set of transducers, respectively tuned to uji2 
and 0^22 ■ The PHC proposal is best described graphically in Figure ^: a second set of resonators, tuned to the second 
quadrupole harmonic a;22 can be placed in an equivalent position in the "southern hemisphere" , and an eleventh 
resonator tuned to the first monopole frequency wio is added at an arbitrary position. It is not difficult to see, by the 
general methods outlined earlier on in this paper, that cross interaction between these three sets of resonators is only 
second order in rj^^^, therefore weak. 

A spherical GW detector with such a set of altogether 1 1 transducers would be a very complete multi-mode multi- 
frequency device with an unprecedented capacity as an individual antenna. Amongst other it would practically enable 
monitoring of coalescing binary chirp signals by means of a rather robust double passage method |38], a prospect 
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FIG. 5. To the left, the pentagonal hexacontahedron shape. Certain faces are marked to indicate resonator positions in a 
specific proposal — see text — as foUows: a square for resonators tuned to the first quadrupole frequency, a triangle for the 
second, and a star for the monopole. On the right we see the (penatgonal) face of the polyhedron. A few details about it: the 
confluence point of the dotted lines at the centre is the tangency point of the inscribed sphere to the PHC; the labeled angles 
have values a = 61.863°, /3 = 87.205°; the angles at the T-vertices are all equal, and their value is 118.1366°, while the angle at 
P is 67.4536°; the ratio of a long edge (e.g. PTi) to a short one (e.g. T1T2) is 1.74985, and the radius of the inscribed sphere 
is twice the long edge of the pentagon, R — 2 PT\ . 



which was considered so far possible only with broadband long baseline laser interferometers [^9 40 1, and is almost 
unthinkable with currently operating cylindrical bars. 



VII. A CALIBRATION SIGNAL: HAMMER STROKE 



This section is a brief digression from the main streamline of the paper. We propose to assess now the system 
response to a particular, but useful, calibration signal: a perpendicular hammer stroke. After a few general consider- 
ations we describe the situation in a TIGA and in a PHC configuration. 

s) in its rhs with that corresponding to a hammer stroke, 



We first go back to equation (3.15 ) and r eplac e 



"external J 



which is easily calculated thanks to ( 3.4c ), ( 3.5c ), and ( |3.16 ); the result is 



"Stroke 



/o 



|A„i(^)r fl(na-no) , 



a = 1, . 



(7.1) 



where no are the spherical coordinates of the hit point on the sphere, and /o = no-io/M. Clearly, the hammer stroke 
excites all of the sphere's vibration eigenmodes, as it has a completly flat spectrum. 

The coupled system resonances are of course those already calculated in section 6, and the same procedures described 
there for a GW excitation can be pursued now to obtain 



(c) (c)* 



P,(n6-no)+0(0) , 



a — \ . . . , J 



(7.2) 
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FIG. 6. Compared line spectrum of a coupled TIGA and a PHC resonator layout in an ideally spherical system. The weakly 
coupled central frequency in the TIGA is drawn dashed. The frequency pair is 5-fold degenerate for this layout, while the two 
outer pairs of the PHC are doubly degenerate each, and the inner pair is non-degenerate. Units in abscissas are rj^^^Q, and 
the central value, labeled 0.0, corresponds to Q. 



when the system is tuned to the n/-th spheroidal harmonic, i.e., Q = u!ni. It is immediately seen from here that the 
system response to this signal when the resonators are tuned to a monopole frequency is given by 



9. 



-4= l^«0(i?)| \ 

V4:nJ ^ 



^ — w„o 



(7.3) 



an expression which holds for all a, and is independent of either the resonator layout or the hit point, wh ich i n 
particular prevents any determination of the latter, as obviously expected. The frequencies uj± are those of (3.12), 
and we find here again a global factor J~^/^, as also expected. 

We consider next the situation when quadrupole tuning is implemented, ft — a;„2 ■ We shall however do so only for 
the TIGA and PHC configurations, as more general considerations are not quite as interesting at this point. 



A. TIGA response to a hammer stroke 



One easily verifies that, for the TI configuration. 



9.(s) = -,r'/2. 



/|^/oK.(i?)|i 



P2{na-no) , TIGA 



(7.4) 



where uj± are the TIGA frequencies ( 6.2Cl| ), and n = uin2- This equation shows that the hitting point position no can 
be easily determined from the readouts 9a(s) — as a matter of fact it is redundantly determined by them. The mode 
channels for this signal yield 



= -rj^y^f, |A„2(i?)| ^ \{s'+coiy' - {s^+uiy' 



^2m(no) , m 



(7.5) 



and they are proportional to the sphere's quadrupole radial oscillation amplitude A„2(^) 5^2m('^o) the hit point, 
no |llo| . A simple numerical simulation is illustrative of the situation, and we present it here. 
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We recall that the inverse Laplace transform of 2 



sin-(w+ cosrJi + 0(?7^/^) 



(7.6) 



which is a beat — a sinusoid of carrier frequency Q amplitude modulated by another sinusoid of much smaller 
frequency of order ry^/^fi. All six resonators' motions are thus identical beats, except that amplitude varies from one 
to the other, and the same applies to all five mode channels. We see this graphically in Figure |^, wherein plots are 
displayed of the resonator readouts and mode channels. 

B. PHC response to a hammer stroke 

The situation is slightly more involved for a PHC distribution, but still easy to address; the system response is 
given by 



m — — 2 



C/i"2™(na)r;„(no) , PHC (7.7) 



and the mode channels by 



(7.8) 



The difference with the TICA is this: the system response qa{t) is a superposition of three different beats, while 
the mode channels are single beats each, but with differing modulation frequencies. This is represented graphically 
in Figure ^, where we see the result of a numerical simulation of the PHC response to a hammer stroke, delivered to 
the solid at the same location as in the TIGA example of the previous subsection. The readouts gait) are somewhat 
fancy time series, whose frequency spectrum shows three pairs of peaks — in fact, the lines in the ideal spectrum of 
Figure |^. The mode channels on the other hand are pure beats, whose spectra consist of the individually separate 
pairs of previous peaks. It might perhaps be said that the PHC gives rise to a sort of "Zeeman splitting" of the 
TICA degenerate frequencies, which can be attributed to an axial symmetry breaking of the isotropic character of 
that resonator distribution; the PHC mode channels naturally resolve the split multiplet into its components. 

This simple example very neatly displays the most relevant features of both TIGA and PHC resonator distributions, 
and enables a quite clear comparison of the merits of either, which does apply also to real GW signals. Despite their 
distinct characteristics, these layouts are perfectly equivalent as regards their ability to sense GWs in an ideally 
noiseless detector. This is so simply because the detector is the sphere, not the transducers, and therefore any motion 
sensing system, not just these two, is equivalent to them if perfectly noiseless monitoring were possible. Things do 
however change in the presence of noise. T. Stevenson has recently addressed the problem of how isotropy in sensitivity 
is affected by noise [|4|, and concludes that while the TICA maintains isotropic sensitivity, the PHC does not: the 
latter is slightly more sensitive in a relatively small solid angle around the resonators' axis, but slightly less for other 
incidence directions. This is a clear disadvantage of the PHC, but one may not forget that a real antenna will not 
be spherically symmetric because it must be suspended, and that breaks that symmetry. As we shall shortly see, the 
PHC is very naturally adapted to that symmetry breaking, and therefore it is probably unwise to outrightly dispose 
of it on the basis of a single noise criterion. At the same time, the PHC system spectrum is richer, and slightly 
broader than the TIGA^s, and this again favours PHC, even in the presence of noise, as the possibility of looking 
for signals at more frequencies reduces the probability of random disturbances, thereby increasing the probability of 
signal detection. 

The above is a qualitative discussion, whose quantitative aspects will not be addressed in any more depth in this 
paper. We still wish to assess the effects of imperfections in a noiseless detector, which are of paramount relevance in 
a real system. 



VIII. SYMMETRY DEFFECTS 



So far we have made the assumption that the sphere is perfectly symmetric, that the resonators are identical, that 
their locations on the sphere's surface are ideally accurate, etc. This is of course unrealistic. So we propose to address 
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now how such departures from ideahty affect the system behaviour. As we shall see, the system is rather robust, in a 
sense to be made precise shortly, against a number of small deffects. 

In order to quantitatively assess ideality failures we shall adopt a philosophy which is naturally suggested by the 
results already obtained in an ideal system. It is as follows. 



As we have seen in previous sections, the solution to the general equations (3.15) must be given as a perturbative 
series expansion in ascending powers of the small quantity 77^/^. This is clearly not a fact related to the system's 
symmetries, and so it will survive symmetry breakings. It is therefore appropriate to parametrise deviations from 
ideality in terms of suitable powers of r]^^^, in order to address them consistently with the order of accuracy of the 
series solution to the equations of motion. An example will better illustrate the situation. 



In a perfectly ideal spherical detector the system frequencies are given by equations (5.4). Now, if a small departure 
from e.g. spherical symmetry is present in the system then we expect that a correspondingly small correction to those 
equations will be required. Which specific correction to the formula will actually happen can be qualitatively assessed 
by a consistency argument: if symmetry deffects are of order 77^/^ then equations (pA) will be significantly alter ed in 



their ry^/^ terms; if on the other hand such deffects are of order 77 or smaller then any modifications to equations (5.4) 
will be swallowed into the 0(0) terms, and the more important r/^^"^ terms will remain unaffected by this symmetry 
failure. We will say in the latter case that the system is robust to that ideality breaking. 

More generally, this argument can be extended to see that the only system deffects standing a chance to have any 
influences on lowest order ideal system behaviour are deffects of order rj^^^ relative to an ideal configuration. Deffects 
of such order are however not necessarily guaranteed to be significant, and a specific analysis is required for each 
specific parameter in order to see whether or not the system response is robust against the considered parameter 
deviations. 

We therefore proceed as follows. Let P be one of the system parameters, e.g. a sphere frequency, or a resonator 
mass or location, etc. Let Pidoai be the numerical value this parameter has in an ideal detector, and let Picai be its 
value in the real case. These two will be assumed to differ by terms of order 77^^^, or 

-Prcal = Pidcal(l+P?7'^') (8.1) 



For a given system, p is readily determined adopting (8.1) as the definition of Preai, once a suitable hypothesis has 



been made as to which is the value of Pidcai- In order for the following procedure to make sensible sens e it is clearly 



required that p be of order 1 or, at least, appreciably larger than ry^/^. Should p thus calculated from {SA) happen 
to be too small, i.e., of order 77^/^ itself or smaller, then the system will be considered robust as regards the affected 
parameter. 

We now apply this criterion to various departures from ideality. 



A. The suspended sphere 



An earth based observatory obviously requires a suspension mechanism for the large sphere. If a nodal point 
suspension is e.g. selected then a diametral bore has to be drilled across the sphere ||2l[]. The most immediate 
consequence of this is that spherical symmetry is broken, what in turn results in degeneracy lifting of the free spectral 
frequencies w^j, which now split up into niultiplcts ujnim 

(m — ~l,...,l). The resonators' frequency cannot therefore 
be matched to the frequency iOn„ig, but at most to one of the LOngiom^S- In this subsection we keep the hypohesis that 
all the resonators are identical — we shall relax it later — , and assume that falls within the span of the multiplet of 
the ^noiom^^- Then we write 



{l+PmV^^^) , m = -lo 



■2) 



We now search fo r th e coupled frequencies, i.e., the roots of equation (3.18). The kernel matrix Kab{s) is ho wever 
no longer given by (4^), due the removed degeneracy of oJni, and we must stick to its general expression ( 3.17 ), or 



0/ I 1 



nlm ' nLm nlm 

Following the steps of section 4 we now need to seek the roots of the equation 



■Xah 



(8.3) 



det 



la 



E 

m=-lr 



(.2+J72)(.2+^2^^^^J 



{nolom) 



E 



(s^ + n^)is^+u;l,J 



(nhn) 
Xah 



(8.4) 
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Since relat es to ujnoigm through equation (^^) we see that the roots of (^^) follow again into either of the two 
categories (|5.2| ), i.e., roots close to ±iV, and root s clo se to ±iu}n im {nl^nolo). We shall exclusively concentrate on 
the former now. Direct substitution of the series (5.2a) into (pM yields the following equation for the coefficient xi- 



det 



1 



lo 

E 



(nalom] 



.5) 



This is a variation of (5.3), to which it reduces when pm — 0, i.e., wh en there is full degeneracy. 

The solutions to (^^) no longer come in symmetric pairs, like (5.4). Rather, there are 2Zo+l+^ of them, with 
a maximum number of 2(2Zo+l) non-identically zero roots if J > 2/o +lp°| . For example, if we choose to select the 
resonators' frequency close to a quadrupole multiplet (/o = 2) then (8.5) has at most 5+ J non-null roots, with a 



ma xim um ten no matter how many resonators in excess of 5 we attach to the sphere. Modes associated to null roots 
of ( ^.5[ ) can be seen to be weakly coupled^ just like in a free sphere, i.e., their amplitudes are smaller than those of the 
strongly coupled ones by factors of order 77^/^. 

In order to assess the reliability of this method we have applied it to see what are its predictions for a real system. 
To this end, data taken with the TIGA prototype at LSU]^ were used to confront with. The TIGA was drilled 
and suspended from its centre, so its first quadrupole frequency split up into a multiplet of five frequencies. Their 
reportedly measured values are 



wi2o = 3249 Hz, wi2i= 3238 Hz, t^ia _i = 3236 Hz , CJ122 = 3224 Hz , tJi2 _2 = 3223 Hz 

All 6 resonators were equal, and had the following characteristic frequency and mass, respectively: 

1 



fl = 3241 Hz , 
Substituting these values into (|8.2| ) it is seen that 



V 



1762.45 



Po = 0.2075 , pi = -0.0777 , p-i = -0.1036 , p2 = -0.4393 , p-2 



-0.4650 



(8.6) 



■7) 



(8.8) 



Equation (8.5) can now be readily solved once the resonator positions are fed into the matrices 



Xab '■ Such 

positions correspond to the pentagonal faces of a truncated icosahedron. Merkowitz ||21[] gives a complete account of 
all the measured system frequencies as resonators are progressively attached to the selected faces, beginning with one 
and ending with six. In Figure ^ we present a graph ical display of the experimentally reported frequencies along with 
those calculated theoretically by solving equation ( |8.5| ). In Table 1 we give the numerical values. As can be seen, 
coincidence between our theoretical predictions and the experimental data is remarkable: the worst error is 0.2%, 
while for the most part it is below 0.1%. This should be taken as very strong evidence that our theoretical model 
is correct, since discrepancies between its predictions and experiment are of order 77, as indeed expected — see ( ^.7] ). 
In addition, it is also reported in reference |2^ that the 11-th, weakly coupled mode of the TIGA (highlighted in 
square brackets in Table ||) has a practically zero ampliutde, again in excelle nt a greement with our general theoretical 
predictions about modes beyond the tenth — see paragraph after equation (8.5). 

This i s an encouraging result which motivated us to try a better fit by estimating the next order corrections, i.e., xi 
of (5.2a). As it turned out, however, matching between theory and experiment did not consistently improve. This is 
not really that surprising, though, as M&J explicitly state |2^ that control of the general experimental conditions in 
which data were obtained had a certain degree of tolerance, and they actually show satisfaction that ^--^1% coincidence 
between theory and measurement is comfortably accomplished. But 1% is two orders of magnitude larger than rj 
— cf. equation (|^) — , so failure to refine our frequency estimates to order 77 is fully consistent with the accuracy of 
available real data. 

We do expect that the analytical procedures developed in this paper will be the appropriate ones to assess and 
theoretically understand the system behaviour as rigor in system parameter control is progressively gained. 



This is a mathematical fact, whose proof is relatively cumbersome, and will be omitted here; we just mention that it has 
its origin in the linear dependence of more than 2Zo-l-l spherical harmonics of order Iq. 

These data are contained in reference and we want to express our gratitude to Stephen Merkowitz for kindly handing 
them to us. 
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TABLE I. Numerical values of measured and theoretically predicted frequencies (in Hz) for the TIGA prototype with varying 
number of resonators. Percent differences are also shown. The calculated values of the tuning and free multiplet frequencies are 
taken by definition equal to the measured ones, and quoted in brackets. In square brackets the frequency of the weakly coupled 
sixth mode in the full, 6 resonator TIGA layout. These data are plotted in Figure H. 



Item 


Measured (Hz) 


Calculated (Hz) 


% difference 


Item 


Measured (Hz) 


Calculated (Hz) 


% difference 


Tuning 


3241 


(3241) 


(0.00) 


4 resonators 


3159 


3155 


-0.12 


Free multiplet 


3223 


(3223) 


(0.00) 




3160 


3156 


-0.11 




3224 


(3224) 


(0.00) 




3168 


3165 


-0.12 




3236 


(3236) 


(0.00) 




3199 


3198 


-0.05 




3238 


(3238) 


(0.00) 




3236 


3236 


0.00 




3249 


(3249) 


(0.00) 




3285 


3286 


0.03 


1 resonator 


3167 


3164 


-0.08 




3310 


3310 


0.00 




3223 


3223 


0.00 




3311 


3311 


0.00 




3236 


3235 


-0.02 




3319 


3319 


0.00 




3238 


3237 


-0.02 


5 resonators 


3152 


3154 


0.08 




3245 


3245 


0.00 




3160 


3156 


-0.14 




3305 


3307 


0.06 




3163 


3162 


-0.03 


2 resonators 


3160 


3156 


-0.13 




3169 


3167 


-0.08 




3177 


3175 


-0.07 




3209 


3208 


-0.02 




3233 


3233 


0.00 




3268 


3271 


0.10 




3236 


3236 


0.00 




3304 


3310 


0.17 




3240 


3240 


0.00 




3310 


3311 


0.03 




3302 


3303 


0.03 




3313 


3316 


0.10 




3311 


3311 


0.00 




3319 


3321 


0.06 


3 resonators 


3160 


3155 


-0.15 


6 resonators 


3151 


3154 


0.11 




3160 


3156 


-0.13 




3156 


3155 


-0.03 




oiyi 


Ql on 


n no 




olDZ 


Ql 


U.UU 




3236 


3235 


-0.02 




3167 


3162 


-0.14 




3236 


3236 


0.00 




3170 


3168 


-0.07 




3297 


3299 


0.08 




[3239] 


[3241] 


[0.06] 




3310 


3311 


0.02 




3302 


3309 


0.23 




3311 


3311 


0.00 




3308 


3310 


0.06 












3312 


3316 


0.12 












3316 


3317 


0.02 












3319 


3322 


0.10 



A final word on a technical issue is in order. Merkowitz and Johnson's equations for the TIGA |l^j22l are identical 
to ours to lowest order in rj. Remarkably, though, their reported theoretical estimates of the system frequencies 
are not quite as accurate as ours. The reason for this is probably the following: in M&J's model these frequencies 
appear within an algebraic system of 5+J linear equations with as many unknowns which has to be solved; in our 
model the algebraic system has only J equations and unknowns, actually equations (3.15). This is a very appreciable 
difference indeed for the range of values of J under consideration. While the roots for the frequencies can be seen 
to mathematically coincide in both approaches, in actual practice these roots are estimated, generally by means of 
computer programmes. It is here that problems most likely arise, for the numerical reliability of an algorithm to solve 
matrix equations normally decreases as the rank of the matrix increases. The significant algebraic simplification of 
our model's equations should therefore be considered one of real practical value. 



1. The suspended PHC 



We briefly consider now which would be the effects of symmetry breaking due to suspension in a PHC resonator 
distribution. The natural suspension axis is the symmetry axis of the resonators, so we shall a ssum e that 5 of them 
are attached to the detector around that axis as described in the paragraphs before equation ( |6.30 ). To be specific, 
we shall speculate with a PHC having the numerical parameters of equations ( |8.6| ) and (^^, a hypothetical but 
reasonable conjecture which will enable comparison with the actual TIGA prototype. 

The results are displayed in Figure It is at once apparent that the structure of the frequency pairs in a drilled 
PHC is quite similar to that of the perfect PHC, while the drilled TIGA pairs are qualitatively different from those 
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of the ideal TIGA. This is not particularly surprising given the symmetries of both layouts, and is an indication that 
axially symmetric distributions may be an interesting alternative to more symmetric distributions in a real system, 
as their geometry naturally adapts to the suspension device, and also require fewer transducers. 



B. Other mismatched parameters 

We finally devote a few words to assess the system sensitivity to small mismatches in resonators' masses, locations 
and frequencies. 



Resonator mass mismatches 



If the masses are slightly non-equal then we can write 

Ma = vM{l + ^lari^^^) , a=l,...,J 



(8.9) 



where 77 can be define d e.g . as the ratio of the average resonator mass to the sphere's mass. It is immediately 
obvious from equation ( |8.9| ) that mass non-uniformities of the resonators only affect our equations in second order, 
since resonator mass non-uniformities result, as we see, in correcctions of order rj^^"^ to r?^/^ itself, which is the very 
parameter o f th e perturbative expansions. The system is thus clearly robust to mismatches in the resonator masses 
of the type ( ^ . 



2. Errors in resonator locations 



The same happens if the locations of the resonators have tolerances relative to a pre-selected distribution. For let 
Ua be a set of resonator locations, for example the TIGA or the PHC positions, and let n'a be the real ones, close to 
the former: 



na (1 + v„ , 



1, 



(8.10) 



The values determine the eigenvalues Ca in equation ( |5.4| ), and also they appear as arguments to the spherical 
harmonics in the system response functions of sections 6 and 7. It follows from (^.10) by continuity arguments that 



(n'a) =n„.(n,) + 0(771/2) 

Ca-Ca+0{r,'/') 



(8.11a) 
(8.11b) 



Inspection of the equations of sections 5-7 shows t hat b oth (a and yim(no) always appear within lowest order 
terms, and hence that corrections to them of the type (8.11) will affect those terms in second order again. We thus 
conclude that the system is also robust to small misalignments of the resonators relative to pre-established positions. 



3. Resonator frequency mistunings 



The resonator frequencies may also differ amongst them, so let 

17, = r!(l + 77^/2) , a=l,. 



i.l2) 



To assess the consequences of this, however, we must go back to equation ( 3.18 ) and see what t he coefficients in 
its series solutions of the type (5.2a) are. The procedure is very similar to that of section VIII A, and will not be 
repeated here; the lowest order coefficient XJ- is seen to satisfy the algebraic equation 



det 



J (nolo) r 
Xac Ocb 



= 



.13) 



which reduces to ( |5.3| ) when all the p's vanish, as expected. This appears to potentially have significant effects on 
our results to lowest order in 77^/^, but a more careful consideration of the facts shows that it is probably unrealistic to 
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think of such large tolerances in resonator manufacturing as implied by equation (8.12) in the first place. In the TIGA 
experiment, for example ||2l|] , an error of order 77^/^ would amount to around 50 Hz of mistuning between resonators, an 
absurd figure by all means. In a full scale sphere (~40 tons, ~3 metres in diameter, ^-^800 Hz fundamental quadrupole 
frequency, ij^ 10~^) the same error would amount to between 5 Hz and 10 Hz in resonator mistunings for the lowest 
frequency. This is probably excessive for a capacitive transducer, but may be realistic for an inductive one. With this 
exception, it is thus more appropriate to consider that resonator mistunings are at least of order 77. If this is the case, 
though, we see once more that the system is quite insensitive to such mistunings. 

Summing up the results of this section, we can say that the resonator system dynamics is quite robust to small 
(of order rj^^"^) changes in its various parameters. The important exception is of course the effect of suspension 
drillings, which do result in significant changes relative to the ideally perfect device, but which can be relatively easily 
calculated. This theoretical picture is fully supported by experiment, as robustness in the parameters here considered 
has been reported in reference [ p2[ . 

IX. CONCLUSIONS 

A spherical GW antenna is a natural multimode device with very rich potential capabilities to detect GWs on 
earth. But such detector is not just a bare sphere, it requires a set of motion sensors to be practically useful. It 
appears that transducers of the resonant type are the best suited ones for an efficient performance of the detector. 
Resonators however significantly interact with the sphere, and they affect in particular its frequency spectrum and 
vibration modes in a specific fashion, which has to be properly understood before reliable conclusions can be drawn 
from the system readout. 

The main objective of this paper has been the construction and development of an elaborate theoretical model to 
describe the joint dynamics of a solid elastic sphere and a set of radial motion resonators attached to its surface at 
arbitrary locations, with the purpose to make predictions of the system characteristics and response, in principle with 
arbitrary mathematical precision. 

We have shown that the solutions to our equations of motion must necessarily be given as an ascending series 
in powers of the small "coupling constant" 77, the ratio of the average resonator mass to the mass of the large 
sphere. The lowest order approximation coresponds to terms of order 77^/^ and, to this order, we recover, and 
widely generalise, other authors' results p^ , |33| , ^ , obtained by them on the basis of certain simplifying assumptions. 
This has in particular enabled us to assess the system response for arbitrary resonator layouts, and to search the 
equations for configurations other than the highly symmetric TJGA, whose properties and/or performance may thus be 
comparatively assessed. This search has led us to make a specific proposal, the PffC, which is based on a pentagonally 
symmetric set of 5 rather than 6 resonators per quadruopole mode sensed. This PHC distribution presents a somewhat 
wider frequency spectrum than the TIGA, and has the interesting property that different wave amplitudes selectively 
couple to different detector modes having different frequencies, so that the antenna's mode channels come at different 
rather than equal frequencies. The PHG philosophy can be extended to make a multifrequency system by using 
resonators tuned to the first two quadrupole harmonics of the sphere and to the first monopole, an altogether 11 
transducer set [p3| . 

The assessment of symmetry failure effects, as well as other parameter departures form ideality, has also interested 
us here. This is seen to receive a particularly clear treatment in our general scheme: the theory transparently shows 
that the system is robust against relative disturbances of order 77 or smaller in any system parameters, and provides 
a systematic procedure to assess larger tolerances — up to order 77^/^. The system is shown to still be robust to 
tolerances of this order in some of its parameters, whilst it is not to others. Included in the latter group is the effect 
of spherical symmetry breaking due to system suspension in the laboratory, which causes degeneracy lifting of the 
sphere's eigenfrequencies, now split up into multiplets. By using our algorithms we have suceeded in numerically 
reproducing the reportedly measured frequencies of the LSU prototype antenna with fully satisfactory precision. The 
experimentally reported robustness of the system to resonator mislocations is also in full agrement with our 
theoretical predictions. 

We take this numerical success as a strong indication that our model is correct. Beyond that, though, we also feel it 
sheds much light into the principles of the functioning of a spherical GW antenna with reson ant transducers, as every 



result ultimately follows from very general and fundamental principles, i.e., equations (2.12). In addition, our master 



equations (3.15) are actually simpler than e.g. Merkowitz and Johnson's p9| , p2| , since they are fewer in number, yet 
they contain more information about the system, as they are accurate to arbitrary order in 77. Lowest order solutions 
to these equations are thus already simpler, and this has an obvious practical value. But higher order corrections can, 
and must, be estimated by suitable analysis of those equations. We have not attempted to comprehensively discuss 
fine structure corrections in this paper, which will likely be necessary as experimental techniques improve. We defer 
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them to near future work. 
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APPENDIX A: 



We give in this Appendix a few important properties of the matrix P;(na-nf,) for arbitrary / and resonator locations 
Ua (a=l,...,J) which are useful for detailed system resonance characterisation. 
We first note the summation formula for spherical harmonics ||2^ 

J- ^ 21 + 1 

J2Yirn{na)Yi,n{nb) = ^—Pi{na-nb), a,6=l,...,J (Al) 

m— — l 

with an obvious change of notation in the arguments to the spherical harmonics, and where Pi is a Legendre 
polynomial: 

To ease the notation we shall use the symbol Vi to mean the entire JxJ matrix P;(na-nf,), and introduce Dirac 
kets \m) for the column J- vectors 



Air 
21 + 1 



m = -l,...J (A3) 



Yim{nj) 

These kets are not normalised; in terms of them equation ( |Al[ ) can be rewritten in the more compact form 

/ 

m— — l 

Equation ( |A4| ) indicates that the rank of the matrix Vi cannot exceed (2/ + 1), as there are only {21 + 1) kets |m). 
So, if J > {21 + 1) then it has at least {J — 21 — 1) identically null eigenvalues — there can be more if some of the ria's 
are parallel, as this causes rows (or columns) of Vi to be repeated. 

We now prove that the non-null eigenvalues of Vi are positive. Clearly, a regular eigenvector, |0), say, of Vi will be 
a linear combination of the kets |m): 

I 

vi 10) = \<t>) , 10) = E 1"^) 

m— — l 

where we h ave called the corresponding eigenvalue, since it is a positive number, as we shall shortly prove. If the 



second (AE) is substituted into the first then it is immediately seen that 

I 



(C' Smm' - (m|m')) 0™' = (A6) 



m' — — l 

which admits non-trivial solutions if and only if 
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det (C'<5w-(m|m')) =0 (A7) 

In other words, are the eigenvalues of the {21 + l)x(2/ + 1) matrix (m|m'), which is positive definite because so 
is the "scalar product" {(f>\(f>')- All of them are therefore strictly positive. 
Finally, since the trace is an invariant property of a matrix, and 

J J 

traceiPi) = Y,Pi{^a-na)=Y.l = J (A8) 

a— 1 a— 1 

we see that the eigenvalues ("^ add up to J: 

traceCT'O-ECa^ E (A9) 



[1] R. Forward, Gen. Rel. and Grav. 2, 149 (1971). 

[2] E. Coccia, G.Pizzella, F.Ronga (eds.). Gravitational Wave Experiments, Proceedings of the First Edoardo Amaldi Confer- 
ence, World Scientific, Singapore, 1995. 

[3] E. Coccia, in M. Francaviglia, G. Longhi, L. Lusanna, and E. Sorace (eds.). General Relativity and Gravitation, Proceedings 
of the GR-14 Conference, World Scientific, Singapore, 1997. 

[4] Ashby N. and Dreitlein J., Fhys. Rev. D12, 336 (1975). 

[5] R.V. Wagoner and H.J. Paik in Experimental Gravitation, Proceedings of the Pavia International Symposium, Acad. Naz. 

dei Lincei (1977). 
[6] There are collaborations in USA, Brazil, Holland, and Italy. 

[7] E. Coccia, V. Fafone, G. Frossati, J. A. Lobo, J. A. Ortega, The hollow sphere as a detector of GWs, submitted to Rhys 
Rev D (1997). 

[8] N.S. Magalhaes et al., Mon. Not. R. Astron. Soc. 274, 670 (1995). 
[9] M. Brunetti, PhD Thesis, Roma (1997). 
[10] J.A. Lobo, Rhys. Rev. D52, 591 (1995). 

[11] Bianchi M., Coccia E., Colacino C.N., Fafone V. and Fucito F. Glass, and Quantum Grav. 13, 2865 (1996). 
[12] Coccia E., Lobo J.A. and Ortega J.A., Rhys Rev D 52, 3735 (1995). 

[13] B.F. Schutz, in Gleiser R.J., Kozameh C.N. and Moreschi O.M., General Relativity and Gravitation, Proceedings of the 

GR-13 Conference, lOP, Bristol, UK (1993). 
[14] P. Rapagnani, Nuovo Gimento 5C, 385 (1982). 
[15] J-P. Richard, Rhys. Rev. Lett. 521, 165 (1984). 
[16] H.J. Paik, in reference above. 
[17] P. Astone et al, Rhys. Rev. D47, 2 (1993). 
[18] Astone P. et al, Europhys. Lett. 16, 231 (1991). 
[19] W. Johnson and S.M. Merkowitz, Rhys Rev Lett 70, 2367 (1993). 
[20] S.M. Merkowitz and W. Johnson, Rhys Rev D 51, 2546 (1995). 
[21] S.M. Merkowitz, PhD Thesis, Louisiana State University (1995). 

[22] S.M. Merkowitz and W. Johnson, to appear in the December-15-1997 issue of Rhys Rev D. 

[23] J.A. Lobo and M.A. Serrano, Europhys. Lett. 35, 253 (1996). 

[24] J.A. Lobo and M.A. Serrano, Class, and Quantum Grav. 14, 1495 (1997). 

[25] Landau L.D. and Lifshitz E.M., Theory of Elasticity, Pergamon Press, Oxford (1970). 

[26] Edmonds A.R., Angular Momentum in Quantum Mechanics, Princeton Univ. Press 1960. 

[27] We use here the better adapted notation a„o and a„2 instead of a„ and &„ , respectively, of reference . Numerical values 

are however identical. 
[28] Brans C. and Dicke R. H., Rhys. Rev. 124, 925 (1961) 
[29] F.G. Tricomi, Integral equations, Interscience Publishers (1957). 

[30] D. Porter and D.S.G. Stirling, Integral equations: a practical treatment, from spectral theory to applications, Cambridge 
Univ. Press (1990). 

[31] C.W. Helstrom, Statistical Theory of Signal Detection, Pergamon Press, Oxford (1968). 
[32] B.H. Bransden and C.J. Joachain, Physics of atoms and molecules, Longman (1984). 



26 



[33] N.S. Magalhacs, O.D. Aguiar, W.W. Johnson and C. Frajuca, GRG 29, 1509 (1997). 
[34] T.R. Stevenson, to appear in Phys. Rev. D, July 15 issue (1997). 

[35] Typically, such formulas contain higher order terms in rj which are calculated on the assumption that only tuned modes 

couple to resonators. They are therefore in error by relative magnitudes of order r?^''^. 
[36] A. Holden, Formes, espa.ce et symetries, CEDIC (1977). 
[37] Z. Har'El, GeometncB Dedicata 47, 57 (1993). 
[38] E. Coccia and V. Fafone, Phys. Lett. A 213, 16 (1996). 
[39] A. Krolak, J.A. Lobo and B.J. Meers, Phys. Rev. D43, 2470 (1991). 
[40] A. Krolak, J.A. Lobo and B.J. Meers, Phys. Rev. D47, 2184 (1993). 



27 



o 

O 



o 
o 



o 
o 
o 



o 

00 



N 

X 

O o 

o e 

1- u 



o 
oo 



CM 



CM 













CM 




! I 



I 1 




T3 
C 
O 
o 



a 



FIG. 7. Simulated response of a TIGA to a hammer stroke: the time series and their respective spectra, both for direct 
resonator readouts and mode channels. Asymmetric widening of spectral lines is due to frequency leakage caused by finite 
integration time. 
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FIG. 9. The frequency spectrum of tfie TIGA distribution as resonators are progressively added from none to 6. Continuous 



lines correspond to measured values, and dashed lines correspond to their rj^^^ theoretical estimates of equation {i.l 
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FIG. 10. A comparative display of the line spectra of a TIGA and a PHC, both in a hypothetically perfect solid shape 
(centre) and in a more realistic, drilled- for-suspension device (bottom and top, respectively). TIGA data are real, see Figure 
bl while PHC data are speculative on the basis of reasonable assumptions — see text. 
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